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We present a simple picture of the gauge factor (GF) enhancement in highly heterogeneous materials 
such as thick-film resistors. We show that when the conducting phase is stiffer than the insulating 
one, the local strains within this latter are enhanced with respect to the averaged macroscopic 
strain. Within a simple model of electron tunneling processes, we show that the enhanced local 
strain leads to values of GF higher than those expected for a homogeneous system. Moreover we 
provide formulas relating the enhancement of GF to the elastic and microstructural characteristics 
of TFRs. 

PACS number(s) 72.20.Fr, 72.80.Tm, 62.20.Dc 
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I. INTRODUCTION 

Thick-film resistors (TFRs) are successfully used as 
piezoresistive sensors due to their high strain-sensitivity 
resistance. Their gauge factors (GFs), conventionally de- 
fined as GF = SR/Re, where SR is the change of resis- 
tance R for an applied strain e, range in fact from GF ~ 2 
up to GF ~ 35 at room temperature!!! Of course, higher 
values of GF lead to a better piezoresistive response of a 
TFR, so that the knowledge of the factors which enhance 
GF is of primary importance in improving piezoresistive 
sensor performances. 

Among these factors, the microstructure of TFRs plays 
certainly an important role in enhancing GF. A typical 
TFR is a composite material in which metallic grains are 
embedded in an insulating glassy matrix, and a positive 
correlation between the mean metallic grain sise And GF 
has been established already since some time.Ercl Fabri- 
cation processes such as the peak firing temperature Tf 
modify the compositional and microstructure properties 
leading to important and complex modifications of GF. 
Typically, GF initially increases with Tf and after having 
reached a-maximum, it shows a more or less pronounced 
decrease.!!! Other important intrinsic factors affecting GF 
are the chemical compositions of both the metallic and 
insulating phases. Finally, in addition to these effects, 
there exists an empirical relationship between (sheet) re- 
sistance R and GF: 



GF cx lni?, 



(1) 



which has been established for room temperatures. Ac- 
cording to this relation therefore the most resistive TFRs 
show also the higher GFs. 

Although these phenomenological properties are useful 
in improving the strain sensitivity in TFRs, the micro- 
scopic physical mechanisms enhancing GF are unclear, 
and the quest for high GF values rests only on empiri- 
cal basis. An even more annoying aspect of this situa- 
tion is that simple tunneling-likc formulas for R are com- 



pletely inadequate to explain the observed GF values. 
For example, at sufficiently high temperatures, transport 
is thought to be governed by tunneling between-adjacent 
metallic grains separated by a mean distance dia 



i?~ R a exp(2d/£), 



(2) 



where £ = h/ y2mV is the localization length and m 
and V are the electron mass and the tunneling barrier 
potential, respectively. The prefactor in Eq.(||) is R = 
<&i?H where i?n = "lirh/e 2 ~ 26 Kf2 is the Hall resistance 
and $ is a dimcnsionlcss geometrical factor. Assuming 
homogeneity, the main dependence on an applied strain e 
comes from the exponent of Eq.(||) through d — ► <i(l + e), 
and the gauge factor reduces to: 



OF = ^ ~ 2d,t 



(3) 



Equations (g) and (g|) are fully consistent with the em- 
pirical relation in Eq.(Q), and R can be rewritten as 
R ~ i?o ex p(GF). However, typical values of the sheet 
resistance are of order 10 — 1000 Kf2, while for example 
for GF = 20 we get R ~ i? o exp20 ~ 10 10 KSV. This 
astronomically large value of R obtained by the simple 
tunneling argument gives an idea of the problem concern- 
ing the physical mechanism of GF in TFRs. 

In this paper we present a theory of GF in TFRs which, 
on one hand, reconciles the observed values of GF with 
those of sheet resistances and, on the other hand, clar- 
ifies the key elements responsible for high gauge factor 
values. Basically, we claim that the high GF values are 
due to pronounced compositional heterogeneity of TFRs, 
so that the local strain ei oc can be substantially different 
from the applied macroscopic strain e. In this situation, 
the mean inter-grain distance is modified by the local 
strain, d — ► d(l + ei oc ), so that Eq.(^) acquires a prefac- 
tor: 



GF ~ (2d/£) 



Qo 



(4) 



The value of ei oc /e depends on the elastic properties of 
the metallic and glassy phases as well as on the mean 
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grain size and inter-grain distance. We show that when 
the insulating phase is less stiff than the conducting one, 
eioc can be substantially larger than e leading therefore 
to an enhancement of GF. In the next section we provide 
analytical formulas relating e\ oc /e to the compositional 
and microstructural characteristics of TFRs by using a 
simplified model of electron hopping and strain distribu- 
tion. 



II. ORIGIN AND BASIC EFFECTS OF STRAIN 
HETEROGENEITY 

Heterogeneity is an intrinsic feature of TFRs. As al- 
ready stated in the introduction, TFRs are two-phases 
materials where metallic particles of sizes ranging be- 
tween ~ 200 and ~ 6000 Aare embedded into an in- 
sulating glassy material. Depending on the fabrication 
processes, a certain amount of conducting material is dis- 
persed in the glass and recent x-ray investigations have 
revealed that much smaller metallic particles (with sizes 
of order ~ 20-70 A) are actually dispersed in the glass.S 
At the present, it is not clear whether the electrons tun- 
nel mainly between the large or the small particles, or if 
possible impurity states in the glass play a positive role 
in the observed conductivities. Several, models of trans- 
port in TFRs have been proposedJ3~II3 having however 
the tunneling process as the common element. Hence, in 
the forthcoming discussion we assume Eq.(||) to capture 
the essential physics at least at high temperatures, and 
interpret d as being a typical tunneling distance ignoring 
for the moment whether it refers to the large or the small 
grain mean separation. 

An important cause of heterogeneity is the chemical 
composition of both the metallic and insulating phases. 
The empirical search for high GF values has selected 
some transition-metal oxides containing heavy elements 
(R11O2, Ir02, Bi2Ru20y are typical examples) to be the 
metallic constituents which present the best TFR perfor- 
mances. An interesting aspect of this selection is that 
such metal-oxides are hard materials characterized by 
high values of bulk-moduli B. For example Ru02, which 
has B ~ 270 GPaJij is among the stiffest jaaaterials and 
Ir02 should have comparable values of B.o Conversely, 
the glassy phase has B values typically of order 40-80 
GPa depending on composition. Hence, the metallic 
phase is much stiffer than the insulating one so that there 
is a pronounced microscopical heterogeneity in the elastic 
properties of TFRs. 

A straightforward consequence of this elastic hetero- 
geneity is that TFRs under an applied macroscopic strain 
e should develop highly variable local strain fields. In 
this case, Eq.(|l|) is inadequate because it has been ob- 
tained by assuming that the microscopic tunneling pro- 
cess, Eq.(^), is affected only by the macroscopic strain. 
Hence, a correct treatment of the GF problem in TFRs 
would require the knowledge of the local strains devel- 
oped in a disordered two-phases material with different 



elastic properties. This is a quite difficult task which 
could be attacked only by extensive numerical calcula- 
tions. However this is outside the scope of this paper, 
instead our aim here it is to show that when a system 
is composed by interconnecting stiff conducting and soft 
insulating phases, the local strain within this latter is en- 
hanced compared to the averaged one, and the tunneling 
process gets a higher strain sensitivity. 

To illustrate in the simplest way the local strain en- 
hancement effect, we consider a one-dimensional model 
defined as a periodic arrangement of N conducting seg- 
ments each of length s separated each other by N insu- 
lating segments of length d, so that the total length of the 
one-dimensional array is L — N(s + d). The stress-strain 
relation is simply <j[x) = 2fj,(x)e(x), where ct(x) and e(x) 
are the stress and strain functions which depend on the 
position x along the array. In the two-phases model here 
considered, the elastic function fiix) assumes the values 
Ho or hi depending whether x lies within the insulat- 
ing or the conducting segments. The equilibrium strain 
distribution is governed by the equation da(x)/dx = 
from which it is found that e{x) is constant within each 
segment assuming the value e(x) = e° (e(x) = e 1 ) for 
x within the insulating (conducting) segments, where 
Hoe = Hie 1 . 

The local values e° and e 1 are related to the macro- 
scopic strain e by the following relation: 

e= T dxe{x) = —- , (5) 

L J Q d+ s 

from which we find that the strain within the insulating 
segments is 

= 1 + s/d 

€ e i + (s/d)(H /HiV {> 
while within the conducting segments is 

c i =c 1 + s/d ho 
e i + {s/d){H /Hi) Mi" 

Equations (Q) and (^) merely illustrate that, if Ho 7^ Mi> 
the strain is concentrated more within the softer phase. 
For example, in the limiting case where the conducting 
segments are perfectly rigid (hi ~ * 00 ) we nn d e 1 = 
and e° = e(l + s/d). 

To investigate the effect of this heterogeneous strain 
distribution on transport, we model the electron tunnel- 
ing between two adjacent conducting segments with a 
tunneling probability proportional to exp(2d/£), so that, 
if we neglect the contribution from the metallic segments, 
the total resistance reduces to: 

R = NRoexp(2d/Q. (8) 

An applied macroscopic strain e modifies the tunneling 
distance d through the local strain, d — > d(l + e ), so that 
the gauge factor GF = SR/eR is 
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GF = (2d/0~. (9) 

By substituting Eq.(g) in Eq.(§), we then obtain 

1 + s/d 



GF = (2d/0 



l + (s/<i)(Mo/Mi)' 



(10) 



The above expression clearly illustrate the enhancement 
effect of heterogeneity driven by both geometrical, s/d, 
and elastic, fj,o/^ii parameters. In fact, for metallic 
phases stiffer than the insulating one, [i\ > fio, GF is 
always larger than 2g?/£, reaching the limiting value 



GF = (2d/0 



Mi 
Mo' 



(11) 



in case of tunneling between large conducting phases 
(s/d 3> 1). Another limiting case of Eq.(|ic|) is given for 
perfectly rigid conducting elements (/ii/fio — > oo) which 
leads to 



GF = (2d/£)(l + s/d), 



(12) 



so that the enhancement factor grows linearly with s/d. 

The simplicity of Eq.(|i"o|) stems from the assumption 
of one-dimensionality for which only one component of 
the stress and strain functions is needed. The analysis 
of a more realistic three dimensional case is of course 
much more complex. In fact, in addition to the complete 
solution of the strain field, one should also calculate the 
actual current path which is given by hopping processes 
along directions also different from the one where the 
field is applied. However, as we show in the following 
section, some approximated analytical results can also 
be obtained for the three-dimensional case by employing 
simplified models for the microstructure of TFRs. 



III. THREE-DIMENSIONAL CASE 

To analyze the piezoresistive response of a TFR let us 
consider a cantilever beam in the x-y plane with its main 
axis lying on the x direction. On top of the cantilever 
is deposited the TFR with its thickness measured in the 
z direction. A bending of the cantilever beam produces 
a strain e in the x direction and, if the thickness of the 
cantilever is sufficiently small, no strain in y. If is the 
Poisson ratio of the cantilever, the strain along the z 
direction is —eu/(l — v). The complete transfer of these 
strain values to the TFR leads to: 



0, e 2 



-ev/{l-v), 



(13) 



where e xx , e yy , and e zz are the macroscopic strain com- 
ponents within the TFR along the x, y, and z directions. 
Within this setup, the piezoresistive response is charac- 
terized by two distinct gauge factors: the longitudinal 
gauge factor, GF^, obtained when the potential differ- 
ence is applied along the x direction, and the transversal 
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FIG. 1. Simplified scheme of the transport pattern within 
the lattice model. The spheres represent the conducting 
phase embedded in the insulating medium (lighter region). 
Each sphere has diameter s and the tunneling processes (in- 
dicated by the solid lines) take place between the shortest 
sphere-to. sphere distance d. 



gauge factor, GFt, when the field is applied along the 
y direction. It is important to distinguish between in- 
trinsic and geometric contributions to the piezoxesistive 
response. In fact, from quite general arguments,Ej it can 
be shown that, under the strain field of Eq.(|T3"|), GF^ 
and GFt satisfy the following relation: 



GF, 



GFt = 2 



G pmtr _ QY T 



intr 



(14) 



where the factor 2 stems from the geometrical deforma- 
tion of the TFR, while GF£ ltr and GF^ tr are intrinsic 
gauge factors which are governed by microscopic trans- 
port properties and are equal for isotropic systems. It is 
however experimentally observed that GF^ — GFt > 2, 
suggesting that the piezoresistive response of TFRs is af- 
fected by a certain amount of anisotropy. 

To describe in a simple way the complex pattern of 
the current path and the local strain distribution, we 
introduce a model in which a cubic lattice of conduct- 
ing spheres of diameter s with near-neighbour distance 
d is immersed in the insulating medium (Fig.l). More- 
over, the two phases have different isotropic elastic prop- 
erties. As already assumed in the previous section, we 
disregard the contribution to transport coming from the 
conducting phase and assume for the moment direct tun- 
neling between two neighbouring spheres. Within this ar- 
rangement of the conducting/insulating phases, we model 
the geometry of the current path and the anisotropy 
by assuming that, in average, charge transport can be 
described by iVn tunneling processes parallel and N± 
tunneling processes orthogonal to the direction of the 
voltage drop.t£l Therefore, if the electric field is applied 
along the x direction, the resulting tunneling resistance 
is R x = N\\Rx + N± (R y + R z )/2, where R x , R y and R z 
are resistances for tunneling along the x, y, and z direc- 
tions, respectively. We make the further approximation 
that R x , R y , and R z have the same average value. When 
the macroscopic strain field of Eq.(|I"3|) is applied, the 
different resistances respond to the corresponding local 
strains: 
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SRx 
R x 



2d 

T 



SRy 



2d o 

p yy 



SRz 
R z 



2d 

T 



(15) 



where e xx — S x d/d, e° y = 5 y d/d and e° z = S z d/d Ac- 
cording to the previous definition and to Eq.(ftq), the 
longitudinal gauge factor reduces to: 



GF 



intr 



SR, 



N\\5R X + N±(6Ry + 5R z )/2 
Ry)/2] 



e xx [N\\R x + N_l(R. 
2d r 

( I \ 

e 



XX \ ^xx ^xx 



(16) 



where \ — N±/(N\\ + N±). By following the same lines, 
the intrinsic transversal gauge factor GF^ tr is estimated 
by using the resistance for a field applied along the y 
direction: R y = N\\Ry + N±(R X + R z )/2. Hence: 



GF^ tr = 



SRy 

^xxRy 



NuSR y + N ± {SR X + 6R z )/2 



e xx [N l{ Ry + N X (R 

2d r 
T 



yy 



(17) 



Note that piezoresistive anisotropy, GF£ ltr > GF^ tr , is 
obtained for x < 2/3. 

At this point, the piezoresistive response can be calcu- 
lated once the values of the local strains are known. How- 
ever, also for the simplified model of spheres on a lattice, 
the general solution of the elastic problem is difficult. 
On the other hand, analytical results can be obtained 
by employing some approximations aimed to simplify as 
much as possible the elastic behavior of such a model. 
In the following we present two different approximation 
schemes. 



A. Rigid spheres 

In our model of transport described above, the current 
path is made of tunneling processes between the shortest 
distances (d) between two neighbouring spheres. Hence, 
the piezoresistive response is governed by variations of d 
along the x, y and z directions due to the applied macro- 
scopic strains of Eq. (|l3|) . These modified inter-sphere 
distances are easily found in the limiting case of perfectly 
rigid spheres since in this case there is not deformation 
of the conducting phase [5 x s — 5 y s = 6 z s = 0). In fact, 
if L x , L y and L z are the linear dimensions of the TFR, 
we have 



SL* 



5 T d 



(18) 



and 



(19) 



and the corresponding expressions for e yy and e zz . Upon 
substitution in Eqs.(nf|,n7h, we find therefore: 



and 



Qpintr 



GF^ tr 



2d 
J 



1 



X2 



2 1 - v 



2d f X l 



2v 



2 1 



(1 + s/d), 



(l + s/d). 



(20) 



(21) 



Although the tortuous character of the current path en- 
ters in the three-dimensional case (through the parameter 
x), equations ( pfj| ) and ( pi] ) have the same enhancement 
factor 1 + s/d as in the one-dimensional case of Eq.([lT 
in the limit fii/fXo — » oo. 



B. s/d > 1 limit 

Another case which can be treated easily is given by 
assuming very short distances d compared to the size of 
the conducting spheres (s/d 3> 1). In this limit and in 
the region where tunneling takes place, the boundaries 
separating two neighbouring spheres can be treated in 
first approximation as planes parallel to each other. If 
we neglect the shear components of the stress and strain 
fields, however expected to be weak in this region, we 
find that the stress and strain fields are locally constant 
within each phase and equal toEJ 



<x \ 


/ Ai + 2/i, Xi 


Xi \ 


f e* 

/ ^XX 


a ly = 


X t X t + 2fj,i 
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1 yy 


<J 


V Xi Xi X 


+ 2 Ml 1 


\ ^zz 



(22) 

where Ai and ^ are the Lame coefficients for the insulat- 
ing (i = 0) and conducting (i = 1) phases. Since at the 
boundaries between the two phases the stresses must be 
equal, we find from Eq.([22|) 




where 



a (3 (3 
= | [3 a (3 

S3 13 a 



(Ai + 2jiti)(Ao + Mo) - AiA 




(23) 



(A + 2 Mo )(A 

Bi ' 
B 



1 - 2 



v\ - ISO 



{l + v{){l-2v Q ) 



(24) 



and 



MoAi — /xiAo 
(A + 2^ )(Ao + Mo) - 
Bi v\ - vq 
Bo (l + i/0(l -2i/ )' 



(25) 
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In the above expressions we have expressed the param- 
eters a and (3 in terms of the bulk modulus Bi and the 
Poisson ratio Vi of the corresponding phase, which are 
linked to the Lame coefficients by the relations: 



3 1 - 



2za 



(26) 



Note that if the spheres are almost touching, the applied 
strain field ( |l5|) is transfered mostly to the conducting 
phase and, if we make the additional assumption that 
within the spheres the strains are almost constant, we 
obtain 



(27) 



which is an approximate but indicative relation. Finally, 
substituting Eqs.(|§-|3) into Eqs.([l6|jl^), we find: 




GF mtr 



2d 



X 2-v f 3 \ 2-v 



v\ - vq 



(l + i/l)(l-2n)) 
for the longitudinal gauge factor and 



-DO 



GF!? tr 



2d 

T 



V 



(l + n)(l-2z/ ) 



1 - 2v 
\-v 

Bo 



IT ( 29 ) 



for the transversal one. 

The main result of Eqs.( [28|]29| ) is that the gauge fac- 
tors are proportional to the ratio of the bulk moduli 
Bi/Bq which, as long as B\ > Bq, determines the en- 
hancement of the piezoresistive response with respect to 
the homogeneous limit B\ = Bq. In Ru02-based TFRs, 
this enhancement factor can in principle be as large as 
Bi/Bq ~ 4—5. Note however that hardening of the glassy 
phase due to diffusion of Ru particles (-Br, u — 220 GPa) 
could reduce Bi/Bq lowering therefore the piezoresistive 
response. 

Another interesting feature is the effect of the Pois- 
son's ratios of the conducting (yi) and insulating (uq) 
phases. In fact, since we consider x < 2/3, when v\ > vq 
the longitudinal (transversal) gauge factor is diminished 
(enhanced) with respect to the limit v\ — vq. For v\ < vq 
this situation is reversed. This is explained by noting that 
Poisson's ratio measures the lateral contraction relative 
to the longitudinal extension for a homogeneous mate- 
rial in a tensile stress. Therefore, in the heterogeneous 
system we consider, the phase with higher Poisson's ratio 
will have a more pronounced transversal contraction than 
the phase with lower value of v. This affects the piezore- 
sistive response as schematically illustrated in Fig. 2, 
where a cut along the x-y plane of our model of TFRs is 
shown. For simplicity, we consider that the only nonzero 





FIG. 2. Pictorial explanation of the effect of different Pois- 
son's ratios of the conducting and insulating (i/q) phases. 
It is assumed that the macroscopic strain is only along the x 
direction. The dashed and solid lines represent the contours 
of the elastically deformed conducting spheres when v\ = v$ 
and v\ > vo, respectively. For an electric field applied along 
the x direction, the gauge factor is diminished while for a field 
along the y direction the transversal response is enhanced. 



macroscopic strain component is e xx — e [i.e., v = 
in Eq.(|l3|)]. The dashed and solid lines represent the 
contour of the spheres deformed by e xx when v\ = vq 
and v\ > vq, respectively. Compared to the v\ = vq 
case, when v\ > vq the inter-sphere distances are re- 
duced in the x direction and enhanced in the y direction. 
Therefore, in the limiting case in which all the tunnel- 
ing processes are aligned with the voltage drop (x = 0) 
we expect that the longitudinal GF is reduced while the 
transversal one is enhanced, in agreement with the trend 
given by Eqs.(^8 29). Of course, finite values of x t en d 
to lower this effect which vanishes in the isotropic limit 
X = 2/3. 



IV. DISCUSSION AND CONCLUSIONS 

The results presented above clearly indicate two essen- 
tial elements which favour high piezoresistive responses 
in TFRs: 

(i) the conducting phase must be stiffer than the glassy 
one (Si > B ), 

(ii) the mean linear size s of the metallic particles must 
be much larger than the typical tunneling distance d. 

Both points (i) and (ii) seem to justify some of the 
trends and properties of TFRs. For example, the use of 
hard oxide-metals as constituent of the metallic phase is 
in agreement with (i), while the observed enhancement of 
the gauge factor with the mean metallic grain sizeou is in 
conformity with (ii). Hence, the problem of how the fab- 
rication processes affect the piezoresistive response could 
be approached by investigating the effects on (i) and (ii) . 
For example, the reduction of GF at high peak firing 
temperature T{ could be due to an enhanced dispersion 
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of the metallic phase into the glass. This would enhance 
the effective bulk modulus of the glass leading to a lower 
ratio Bi/B . 

In conclusion, we have proposed a theory of piezoresis- 
tivity in TFRs based on elastic heterogeneity. Within a 
simplified model of TFRs, we have shown that when the 
conducting phase is stiffer than the insulating one, the 
local strains within this latter are enhanced with respect 
to the averaged macroscopic strains. This enhancement 
leads to higher piezoresistive responses leading to higher 
GF values compared to those expected for a homogeneous 
system. 
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we find therefore Bi ~ B 2 . This estimate holds true also 
when more refined relations between 9d and elastic con- 
stants are used [H. Siethoff and K. Ahlborn, J. Appl. Phys. 
79, 2968 (1996)]. 

B. Morten, L. Pirozzi, M. Prudenziati, and A. Taroni, J. 
Phys. D: Appl. Phys. 12, L51 (1979). 

16 C. Grimaldi, P. Ryser, and S. Strassler, J. Appl. Phys. 
(2000); in print. 

17 J. R. Barber, Elasticity (Kluwer, Dordrecht, 1992). 



1 M. Prudenziati, Handbook of Sensors and Actuators (Else- 
vier, Amsterdam, 1994), p. 189. 

2 P. F. Carcia, A. Suna, and W. D. Childers, J. Appl. Phys. 
54, 6002 (1983). 

3 M. Prudenziati, B. Morten, F. Cillone, and G, Ruffi, Sen- 
sors and Actuators 19, 401 (1989). 

4 M. Tamborin, S. Piccinini, M. Prudenziati, and B. Morten, 
Sensors and Actuators A 58, 159 (1997). 

5 Here we disregard additional temperature dependences 
given by thermal expansion, intra-grain electron-phonon 
scattering, activation energies, etc.. These effects can 
be added, however they are not important for order-of- 
magnitude arguments presented here. 

6 C. Meneghini, S. Mobilio, F. Pivetti, I. Selmi, M. Pruden- 
ziati, and B. Morten, J. Appl. Phys. 86, 3590 (1999). 

7 G. E. Pike and C. H. Seager, J. Appl. Phys. 48, 5152 
(1977). 

8 P. Sheng and J. Flatter, Phys. Rev. B 27, 2583 (1983). 

9 M. Prudenziati, A. Rizzi, P. Davoli, and A. Mattei, Nuovo 
Cimento 2, 697 (1983). 

10 W. Schoepe, Physica B 165-166, 299 (1990). 

11 J. Roman, V. Pavlik, K. Flachbart, C. J. Adkins, and J. 
Leib, J. Low Temp. Phys. 108, 373 (1997). 

12 K. Flachbart, V. Pavlik, N. Tomasovicova, C. J. Adkins, 
M. Somora, J. Leib, and G. Eska. Phys. Stat. Sol. B 205, 
399 (1998). 

13 R. M. Hazen and L. W. Finger, J. Phys. Chem. Solids 42, 
143 (1981). 

14 We are not aware of reported measurements of B in RO2. 
However, we can estimate t he bulk modulus of RO2 by 
using the relation 60 — C\J aB/M, where 80 is the De- 
bye temperature, a the lattice constant and M the molec- 
ular mass (C is a constant). Among the oxides, LO2 is 
the compound most similar to RUO2 concerning both the 
electronic and structural properties. In normal conditions, 



G 



